The phase diagram of QCD in a BAG+HRG based equation of state: appearance of a 

pseudo-critical point 
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Mapping the QCD phase boundary and locating critical end point still remains as an open problem 
in strong interaction physics. Predictions about the co-ordinates of the critical point in the {T, hb) 
plane, from different QCD motivated models show a wide variation. Lattice QCD calculations are 
also available, that give an estimation of the critical point for chiral phase transition, where the 
transition changes its nature from rapid cross over to first order transition. Recently co-ordinates of 
the critical point for deconfinement phase transition are claimed to be found as an endpoint of the 
first order phase transition line, in Bag model scenario. In the present paper we have shown that 
Bag model gives a complete first order phase transition line in the (T, /is) plane, and one can not 
have any point where the transition changes its nature. 
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In recent years there is much interest in QCD phase 
diagram and in particular about the location of QCD 
critical end point. That a critical end point exists in the 
QCD phase diagram can be argued as follows: finite tem- 
perature lattice simulation of baryon free {ps = 0) QCD 
indicate a cross-over transition at « 170 ±20 MeV 3, 
1 01 ; • On the other hand there are indications that 
at zero temperature but at finite baryon density, QCD 
has a 1st order transition Q, 0, Q- Only with the 
existence of a critical end point where the 1st order tran- 
sition ends, the cross over transition at (T,fiB — 0) can 
be reconciled with a 1st order transition at (T = 0,/is). 
Unfortunately, lattice simulations at finite baryon den- 
sity is plugged with the well known fermion sign problem 
and at the present state of art, location of QCD phase 
boundary or the critical end point is beyond the lattice 
simulations. Several methods have been devised to sur- 
pass the fermion problem, e.g Taylor expansion Q, [3, 
re- weighting [Tol. [ll| . [l2j . imaginary chemical poten- 
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However, it is uncertain whether the 



methods mentioned above are valid at large baryon den- 
sity. QCD phase diagram has been studied in QCD in- 
spired models, e.g. instanton models, linear sigma model, 
Nambu-Jona-Lasinio (NJL) or Polyakov loop extended 
Nambu-Jona-Lasinio models [1], [3, H, [ll, [ill- 
In QCD inspired model, in general quark degrees of free- 
dom are integrated out. The model results depend on the 
cut off parameter. Depending on the model parameters, 
QCD phase diagram and location of critical end point 
vary widely. 

In a recent work, Singh et. al. [l^], have obtained 
the QCD phase boundary for de-confinement transi- 
tion by using Gibb's equilibrium criterion for first or- 
der phase transition between QGP and hadronic phase. 
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The hadronic phase is modeled by a gas of interacting 
hadrons, where the geometrical size of the baryons is ex- 
plicitly incorporated as the excluded volume correction, 
in a thermodynamically consistent manner. All baryons, 
mesons and their resonances having masses up to 2 GeV 
are included in calculation and strangeness conservation 
is taken into account by equating the net strangeness to 
zero. Mesons are considered as point-like particles, where 
as an equal volume Vex = is assumed for each type of 
baryon with a hard-core radius = 0.6 — 0.8 fm. In ad- 
dition, full quantum statistics is included in the partition 
function of the grand canonical ensemble which helps to 
navigate the region of the phase diagram with low T and 
high /is . Thermodynamic functions of weakly interacting 
quark matter is obtained by a simple Bag model equation 
of state (EOS) with perturbative corrections of the order 

3 /2 • 

of tts in strong interaction coupling constant a^. Hav- 
ing modeled the two phases, the QCD phase boundary is 
obtained advocating Gibb's equilibrium criterion for first 
order phase transition. Gibb's criterion demands that at 
the transition point the pressure of the hadronic state 
Ph{T, 11b) and that of the qgp state Pq{T, iib) becomes 
equal. 

Starting from a low but non-zero value of T and large 
value of /is , gradually they have moved towards large T 
and small /is. At a given T, a corresponding value of 
/is is found out at which the pressure equality holds and 
beyond which QGP pressure dominates. It has been ob- 
served that the line of co-existence between two phases 
ends at a point which is being interpreted as the end 
point of the 1st order transition line, or the critical end 
point. Thus the precise co-ordinates of the QCD critical 
point has been claimed to be estimated for deconfinement 
phase transition in the /is — T plane. The critical val- 
ues of temperature and chemical potential are found to 
be {Tc = l60MeV,^ic = 156MeV) for hard-core radius 

= 0.6/to and Bag constant S^/^ — 216 MeV. For 
Tc = 0.8/m and B^/* = 200 MeV, the point is shifted 
down in temperature to (Tc = lAQMeV, fic = 156MeV). 
Variation in QCD scale parameterA is reported to give 
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insignificant change in the location of the phase bound- 
ary. Considering that earher investigation [2^ , [s^l , dH , 
[22], of phase diagram in terms of Bag model EOS 
and hadronic resonance gas fails to detect the QCD crit- 
ical point, we have studied the model in detail. We find 
that the simple interpretation of minimum of the chemi- 
cal potential below which co-existence line ceases to exist 
as the critical end point is misleading and appears due 
to lack of sufficient numerical precision. 

In the present work we have employed the same Bag 
model EOS as used in [2J| to compute the thermody- 
namic parameters of the partonic phase. In this model, 
QGP is assumed to consist of massless quarks (u,d), their 
antiquarks and gluons only and the pressure in the QGP 
phase takes the form 
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where T dependence of as can be given as [15]: 



as 
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ln{ 



0.089/i| + 15.622r2 



A2 



(2) 



Here we have used B^/^ = 200 MeV and A = 100 MeV 
in our calculation. 

For simulating the hadronic phase we considered a 
hadron resonance gas model (HRG) which includes all 
hadrons and their resonances having masses up to 2.5 
GeV. Excluded volume correction has also been incorpo- 
rated following the prescription by Cleymans and Suho- 
nen (33 |. In this approach, for a given eigen volume Vex, 
the excluded volume corrected pressure, is obtained as, 



Pm^.I^b) (3) 



Here Pgf^j^j-^ [T,^b) is the pressure of the corresponding 
the ideal gas having point like baryons (mesons) and can 
be computed as 



1 



,(^m2+p2_^,„,)/T 



± 1 



(4) 



where the sum is over all the baryonic (mesonic) species 
and their resonances included in the calculation. 

Though not thermodynamically consistent, but the re- 
sulting reduced pressure obtained by this approach is in 



FIG. 1: (color online) The estimated QCD phase boundary 
for first order deconfinement phase transition, with B^''* = 
200 MeV and Vex ~ 1 fm, obtained through both T-scan 
and /^s-scan. The step length for search is kept same for 
both cases (see text). 



close agreement with the approach employed by Singh 
et al (35|, which preserves thermodynamical consistency. 
Having computed the pressure of the two phases, we then 
obtain the phase diagram, in /xs — T plane, by searching 
for zeros of the equation. 



AP{T,^^B) = PQGp{T,fiB)-Pr'{T,fiB) = (5) 

The locus of the zeros is identified as the phase bound- 
ary in hb — T plane along which QGP and the hadronic 
phase can co-exist. The numerical zero-search can be 
employed in two ways. We can fix a temperature Tq and 
vary fj, b over the specified range in steps of Sfi b and cal- 
culate AP(ro 7 A's) a.nd AP{To, ^B + S^b)- If they are of 
opposite sign then there must a value of /is = /i^, such 
that AP{To, = 0, and which can be calculated using 
straight line interpolation. Then (TqiMb) denotes a first 
order phase transition point in the ^b — T plane, pro- 
vided AP(To, /^s (5/^5 ) > AP(To,/xb). We can then go 
to another temperature Tq+ST and get the corresponding 
value of fjf^ . In this way we can trace the first order phase 
transition line in fiB — T plane. Let us name this method 
as ^B-scan. This is identical to the method employed by 
Singh et. al. to construct the QCD phase digram. On 
the contrary one can in principle fix a baryon chemical 
potential /i^ and scan the temperature axis identically 
to obtain the corresponding transition temperatures. We 
can call this approach as T-scan. 

We have employed both methods namely T-scan and 
/i^-scan to draw the phase diagram in the hb —T plane. 
For both cases we have fixed the step length of scan as 
ST — SfiB = 5 MeV. The resulting phase boundaries are 
shown in Fig. [T] 

The first order phase transition line obtained from /is 
scan ends at some non-zero value of /is , the so called crit- 
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FIG. 2: The QCD phase boundary obtained from /is-scan, for gradually decreasing values of 6T. As 5T decreases, the line 
gradually approaches towards temperature axis. For all the cases we have assumed B^^* = 200 MeV and Vex = 1 fm. 



ical point of the de-confinement phase transition. How- 
ever T-scan with same step length yields a phase bound- 
ary that reaches up to /i^ = 0. Hence the end point of 
the first order phase transition line as obtained through 
fiB scan can no longer be interpreted as the true criti- 
cal point of de-confinement phase transition. Rather it 
is a pseudo-critical point, which appears due to insuffi- 
cient numerical resolution of the method employed for the 
construction of the phase diagram. The absence of 1st 
order phase transition points (zeros of AP{T, hb)) be- 
yond some minimum of /i^ can be understood by looking 
at the structure of the phase diagram obtained through 
T-scan. As the line approaches towards the T axis it 
becomes more and more flat in T. Hence the zeros of 
AP becomes more and more closer in T and lie with in 
the search length for a step length of ST — 5 Mev. If one 
gradually decreases the step length, more and more zeros 
start to appear in the low fj.B region and corresponding 
curve eventually moves towards the /is = axis. This 
becomes evident from Fig. [2l where we have plotted the 
phase transition line for gradually decreasing values of 



6T. Thus the appearance of the end point of the first 
order phase transition line in /is — T plane, is not phys- 
ical. This is rather a pseudo-critical point, coming into 
the picture due to lack of numerical detectability. 

Indeed, one can argue that for bag model type equation 
of state, it is possible to tune the bag pressure such that 
Eql5] is satisfed for zero baryons density. In Fig. [31 for 
Bi/4=(50-300) MeV, we have plotted the pressure dif- 
ference AP(T) = Pqgp{T,hb = 0) - Ph{T,ijlb = 0) 
as a fuction of temeprature. One observes that for 
b^MeV < B^/^ < 270 MeV, a co-existence phase ex- 
ists. At low temperatures, QGP pressure is smaller than 
the hadron pressure and AP{T) is negative. At high 
temperature QGP pressure dominates over the hadronic 
pressure resulting a positive AP{T). However, the tran- 
sition temperaure shifts to higher T with increasing B 
and there by indicating a Bag pressure dependent first 
order de-confinement phase transition at zero baryon 
chemical potential. At sufficiently high (low) values of B 
{B^/-^ = 300(50) AfeF), the QGP pressure is always lower 
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FIG. 3: Variation of AP with temperature at zero chemical 
potential for different values of Bag constant. 



FIG. 4: Variation of the Helmholtz free energy density with 
temperature at zero chemical potential for both hadronic 
phase and QGP phase. 



(higher) than the hadronic pressure (at fiB = 0) , and one 
can not get any first order phase transition. But over a 
reasonable range of Bag constants, one can always get a 
solution of AP(T, 0) = and hence a first order hadron 
to QGP phase transition at /zb = 0. At very high values 
of temperature AP{T) is seen to decrease with T. This 
can be attributed to the larger degrees of freedom as- 
sociated with the hadronic phase compared to the QGP 
because of the exponential growth of the hadrons and res- 
onances at very high T, resulting in a higher pressure in 
HG than QGP.This indicates a reversal phase transition 
from QGP to hadron gas at a still higher temperature 
[ssl ]. But once the system goes over to the QGP phase, 
due to rise in temperature, it should continue to stay 
in that phase owing to the asymptotic freedom of QCD. 
Moreover, it is expected that the hadronic interactions 
are become significant when hadrons are closely packed 
in a hot and dense hadron gas. This anomalous behavior 
arises owing to the treatment of the hadronic phase as 
an ideal gas of non-interacting point-like hadrons. As a 
result of this assumption, the thermal production of an 
arbitrarily large number of hadrons in a given volume at 
very high T (or muB) is possible and eventually leads 
to infinitely large energy densities and pressure. In fact 
a simple remedy to this problem is the inclusion of fi- 
nite, proper volume for each hadron, which leads to a 
hard-core repulsion among themselves at very high tem- 
perature and/or density and thereby limiting the number 
of hadrons in the system so that its volume is completely 
filled with particles. This is the so called excluded ex- 
cluded volume correction, where the repulsive force is 
being incorporated by assigning a geometrical hard-core 
volume to each hadron. This finally leads to the reduc- 
tion of the effective pressure and energy density of the 
hadron gas particularly at high temperature and/or den- 
sity. We do have incorporated the excluded volume effect 
in our calculation following the model by Cleymans and 



Suhonen 3J]. In this model the the volume correction 



term is proportional to the net baryon density and thus 
becomes in-effective at muB = 0. Hence the hadron gas 
behaves like an ideal gas and shows the indications of the 
reverse phase transition at very high temperature. 

That our observation of the absence of any true criti- 
cal point in HRG+Bag model based EOS can be justified 
further through following arguments. The critical point 
is the end point of the first order phase transition line 
and at the critical point transition is believed to be sec- 
ond order. A first order phase transition is always associ- 
ated with non-zero finite latent heat, which in fact mea- 
sures the discontinuity in entropy at the transition point. 
For a second order transition, entropy changes smoothly 
along the transition point and hence latent heat is zero. 
In Fig. m we have shown the temperature variation of 
Helmholtz free energy density (/ = e — Ts) at = 0. 
Both fhad and Jqgp decreases as temperature increases. 
Since the system always follows the path of minimum 
free energy, hence the system will undergo a transition 
from hadron phase to partonic phase at T — 129.6 MeV. 
Since s = -■§p, hence a kink in f(T) at T = U^.QMeV 
would indicate a discontinuity in entropy density thereby 
ensuring the zero chemical potential, temperature driven 
transition, to be of first order. 

To conclude, we would like to clarify that we are not 
doubting the existence of the critical point in the QCD 
phase diagram. The ab-initio Lattice QCD calculations 
have indeed proved that such a point does exist. Though 
its precise location in the T — fiB is a matter of ongoing 
debate as predictions from different groups vary wildly. 
Our only motivation is to prove that the explicit con- 
struction of the first order phase separation boundary 
through comparison of a HRG EOS with Bag model EOS, 
can not give any estimation of the critical end point for 
de-confinement phase transition. 
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